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Abstract 

The vacuum of Standard Model above electroweak transition can be regarded as a system with 
marginal Fermi points (MFP) [HE], which is not stable and may decay into three topologically 
distinct universality classes, i.e., vacua with gap, fermi surfaces, or Weyl points. We expect that, 
condensed matters with emerging relativistic quantum field theory in the low energy corner, may 
be relevant to such physics as well. Here we show, based on first-principles calculations, that 
crystalline ^sBi (A=Na, K, Rb) are three dimensional (3D) Dirac semimetals with MFP protected 
by crystal symmetry. They possess non-trivial Fermi arcs on the surfaces, and can be driven 
into various topologically distinct phases by explicit breaking of symmetries. Giant diamagnetism, 
linear quantum magnetoresistance, and quantum spin-Hall eff^ect will be expected for such novel 
compounds. 
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The topological consideration of effective relativistic quantum field theory in 3D mo- 
mentum space allow us to classify quantum vacua (into three distinct classes) and simulate 
high-energy physics and cosmology by probing condensed matters P, 12]. The topological 
classes with fermi surface or gap are well known in condensed matters as metals or insula- 
tors, and the fully gapped systems have been further classified recently into topologically 
trivial and non-trivial insulators [3H5]. The class of compounds with Weyl fermi points are 
rare, however, the A-phase of ^He [6] and recent proposals [THTO] suggested various possi- 
bilities. Given those known realizations of distinct topological states, it is still challenging 
to have a well controlled example near phase boundary such that various topological phase 
transitions can be studied systematically within one system. This challenge becomes further 
relevant because the vacuum of Standard Model is regarded to be at the phase boundary 
with MFP, where two Weyl points with opposite chirality (or topological charge) overlap, 
forming massless Dirac fermionic gas [U |2]. The Weyl points with opposite chirality are 
stable topological objects only if they are separated. If they meet, their topological charges 
may cancel each other and open a gap. In principle, we may obtain accidental MFP by 
fine tuning of chemical composition or spin-orbit coupling (SOC) strength pTHT^ . Unfortu- 
nately, such realizations are too fragile and hard to be controlled. Here we show that, in the 
presence of crystal symmetry, MFP can be protected and stabilized |l5j, and this scenario 
can be achieved in a simple stoichiometry compound A^Bi (y4=Na, K, Rb), which supports 
3D massless gas of Dirac fermions, different with graphene and Weyl semimetals. 

Among the alkali pnictides A^B (y4=Alkali metal, B=As, Sb or Bi), A^Sh are well known 
for its application as photocathode materials [16], but the physical properties of A^Bi are 
not widely studied [171 . Both LiaBi and CsaBi crystallize in cubic Fm3m structure, while 
NasBi, KsBi and RbsBi are in hexagonal PGs/mmc phase (or D^f^, shown in Fig.l) [T8] . 
which are our main focus here by taking NaaBi as an example. In this structure [19], there 
are two nonequivalent Na sites (Na(l) and Na(2)). Na(l) and Bi form simple honeycomb 
lattice layers which stack along the c-axis, while Na(2) atoms are inserted between the 
layers, making connection with Bi atoms. From the ionic picture, due to the closed-shell 
configuration where the number of valence electrons (3xNa-s^+Bi-p'^) is equal to 6, we may 
expect semiconducting nature of these compounds, similar to NasSb [20j. However, they are 
in fact different. 

The calculated electronic structures (see Methods) shown in Fig. 2 suggest that the 
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FIG. 1: (Color online) (a) Crystal structure of NasBi with PQ^/mmc symmetry. Na(l) is at 2b 
position ±(0,0,^), and Bi is at 2c position ib(|,|,^). They form honeycomb lattice layers. Na(2) 
are at 4/ position ±(g, and ib(|,|,2+?i) with ti=0.583, threading Bi along c axis, (b) Brillouin 
Zone of bulk and the projected surface Brillouin Zones of (001) and (010) plane. 



valence and conduction bands are dominated by Bi-6p and Na-3s states. Very close to the 
Fermi level, the top valence band is mostly from 'Qi-Qpx y states, while the conduction band 
with very strong dispersion is mostly from Na(l)-3s states. All these pictures are similar 
to that of NasSb [20], but with the key difference: at the F point, the Na-3s band is lower 
than 'Qi-Qpx,y by about 0.3eV, and it is further enhanced to be 0.7eV in the presence of 
SOC, resulting in a metal with inverted band structure, rather than the normal narrow gap 
semiconductor like NaaSb [20]. The band inversion is mostly due to the heavier Bi, which has 
higher 6p-states and larger SOC compared to Sb. Considering the possible underestimation 
of band gap by GGA, the band inversion can be further confirmed by the following evidences: 
(1) calculation using hybrid functional (HSE) [21] gives band inversion around 0.5 eV, still 
reasonably strong; (2) the earlier calculations for normal semiconductor KaSb [22] suggest 
that its experimental gap can be reasonably reproduced by LDA. 

Having the inverted band structure, however, NasBi is not gapped, different from topolog- 
ical insulators like Bi2Te3 and Bi2Se3 [23] • It is a semi-metal with two nodes (band-crossings) 
exactly at fermi level (Fig. 2). In other words, its fermi surface consists of two isolated fermi 
points, which are located at (0, 0, A;^^±0.26x^) along the F-A line. Since both time re- 
versal and inversion symmetries are present, there are four-fold degeneracy at each fermi 
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FIG. 2: (Color online) The calculated electronic structures of NaaBi. (a) The total and partial 
density of states, (b) and (c) are the band structures without and with spin-orbit coupling, re- 
spectively. The red circles indicate the projection to the Na-3s states. The orbital characters of 
wave- functions at F point are labeled in the inset (see Method for details). 

point, around which the band dispersions can be linearized, resulting in a 3D massless Dirac 
semimetal. It is different from graphene not only in dimensionality, but also in its robust- 
ness, because the fermi points here survive in the presence of SOC. This fact also makes 
difference with other proposals [2H ■ 

The 4x4 Dirac fermion here is massless because the two bands which cross each other 
along the F-A line belong to different irreducible representations under three-fold rotational 
symmetry. Breaking of this symmetry will make the system insulating. For example, 1% 
compression along y will open up a gap ^5.6meV. This insulating state, however, is topo- 
logically non-trivial with ^2=1 OH] due to the inverted band structure around the F point. 
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FIG. 3: (Color online) The projected surface states and their fermi surfaces of NasBi. (a) and (b) 
are the surface states for (001) and (010) surfaces, respectively, (c) The fermi surfaces (fermi arcs) 
and its spin texture (in-plane component) for the (010) surface states, (d) The fermi arcs of (010) 
surface obtained from fitted effective hamiltonian with additional exchange field hi=6 meV (see 
Method for details). The discontinuity around the singular fermi points becomes now obvious. 

This fact makes NasBi unique, because both bulk fermi points and non-trivial surface states 
(a single pair) should coexist (see Fig. 3) as long as the crystal symmetry stands. Further- 
more, the surface states are different from that of topological insulators [23], in the sense 
that its fermi surface has fermi arc structures. As show in Fig. 3(b) for the (010) surface of 
stoichiometry NasBi, although the entire fermi surface is closed, its derivative and fermi ve- 
locity are ill defined at the two singular points (corresponding to the projection of bulk fermi 
points to the surface). The spin texture of surface states has helical structure (also similar 
to topological insulators), but the magnitude of spin vanishes at the singular points. This 
kind of fermi surfaces has never been found before, and it can be understood following the 
discussions for Weyl semimetal |71IH]. If we split the 4x4 Dirac point into two separated 2x2 
Weyl points in momentum space by breaking time reversal or inversion symmetry P, HU] , 
the fermi surface of surface states will also split into open segments which are fermi arcs 
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discussed in Weyl semimetal (as shown in Fig. 3d) [H |8]. All these characters in contrast 
to conventional metals and topological insulators should be experimentally measurable by 
nowadays ARPES technique. Our further calculations for KaBi and RbaBi suggest that they 
have qualitatively the same pictures as NasBi. 



Effective Hamiltonian 



The low energy physics here can be understood from effective Hamiltonian. At the F 
point, the low energy eigen states can be approximated as the linear combination of Na-3s 
and Bi-6p atomic orbitals. Among them, we focus on the four states |S't,+^), 

2 2 

I'S'i , — \Ps , — i), which capture the band inversion nature (see methods, here the super- 

2 2 

script ± indicate the parity of each state). We can then construct an effective 4x4 k ■ p 
Hamiltonian using these four states as basis. Considering the time reversal, inversion, and 



Dqi^ symmetries, the leading order Hamiltonian around F reads: 



Hr{k) = eo(k) + 



M(k) Ak+ B*{k) 

Ak^ -M(k) B*{k) 

B{k) M(k) -Ak^ 

y B{k) -Ak+ -M(k) J 

where eo(k) = Co + Cik^, + ^2(^2 + k^), k± = k^ ± iky, 5(k) = B^k.kl, and M(k) = 
Mq — Mikl — M2{k1 + fcy) with parameters Mq, Mi, M2 < to ensure band inversion. Please 
note the leading order term of off-diagonal elements B(k) has to take the form of k^k^ under 
three fold- rotational symmetry and the opposite parity of IS") and \P) states. Evaluating 



the eigen values -E'(k) = eo(k) ± J M(kY + A^kj^k- + |i?(k)p, we get two gapless solutions 



at k'==(0, kl = ± 



which are the MFP discussed above. There is a pair of such 



points, which are separated along the T — A line. 

If we only concentrate to the neighborhood of each crossing points k'^, and neglect the 
high order terms (i.e., -B(k) ^ 0), the linearized Hamiltonian is nothing but 3D massless 
Dirac fermions. The block diagonal form allows us to decouple the 4x4 matrix into two 2x2 
matrices, which are Weyl fermions with degenerate energy but opposite chirality [THTO]. The 
violation of three-fold rotational symmetry, however, will introduce a linear leading order 
term of -B(k), i.e, B{k)=Bikz. In such case, two Weyl fermions will be coupled together, 
resulting in massive Dirac fermions with gap, similar to the case of Bi2Se3 or Bi2Te3 




FIG. 4: (Color online) Phase diagrams of NasBi with mass term m=0 (left panels) and m=5.6meV 
(right panels). The high order term of -B(k) is neglected for the case of m=5.6meV. (a) and (b) are 
phase diagrams, and (c)-(j) show band dispersions, corresponding fermi surfaces and its topological 
charges for some characteristic phases (with hi fixed to be 10.0 meV). Only the neighborhood 
around one of the MFP is show with the k defined as k — k^. The | ± |), I i |) are abbreviations 
(i.e. Jz values) for the four basis, which are used to indicate the main component of the wave- 
functions for the states away from band crossings, (see text part for details). 

Nevertheless, as long as the three-fold rotational symmetry survives, the MFP here should 
be stable and protected. 
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Phase diagram and topological phase transitions 



NasBi with MFP is just located at phase boundary and may be driven into various 
topologically distinct states by exphcit breaking of symmetries [H, For simphcity of 
illustration, here we focus on the effects of exchange splitting, assuming it can be induced 
by magnetic doping. Other symmetry-breaking terms (such as inversion, mirror or two 
fold rotational symmetries of the crystal) may play the same roles and can be analyzed 
analogously. Since the IS") and \P) are different orbitals (or pseudo spins), we may in 
general separate any exchange splitting into orbital dependent and independent parts as 
Hexi = hiaz®Tz and Hf,x2 = h2(Jz^I, where hi and /i2 are field strength (along z direction), 
and a and r are pauli matrices describing spin and pseudo-spin respectively. The total 
Hamiltonian is given as if = Hr + Hexi + Hex2, and the resulting phase diagram is shown 
in Fig.4. 

If the three fold rotational symmetry of crystal is preserved (i.e. B{\s)=B^klk\^ left 
panels of Fig.4), starting from the MFP state (/ii=0, /i2=0), the state with Weyl points will 
be introduced by hi (horizontal axis), because such exchange field will split the Dirac point 
into two separated Weyl points in momentum space. If the /12 is further introduced, however, 
the two Weyl points will separate energetically, and a system with fermi surfaces (FS) will 
be obtained. On the other hand, if the three-fold rotational symmetry is broken, a mass 
term m ^ Bik'^ will be induced, and the high order term of B{k.) can be then neglected. 
In such case, a topologically non-trivial insulating phase will appear (right panels of Fig.4), 
and the Weyl semimetal phase can be driven only if hi is larger than the mass term m. 

The FS states can be further classified according to their enclosed topological charge 
(or Chern number Cps), which is defined as the net flux of the Berry phase gauge field 
penetrating the fermi surface. 



where the integrand is the Berry curvature, A(k)=— Vkl^^k) is the adiabatic Berry 
connection for the states \uk) at fermi level, and dS points from low to high energy. For the 
case m=5.6meV (right panels), the two distinct FS states, trivial {Cps = 0) and non-trivial 
{Cfs = =tl)) are separated by the line defined as hi — h\=m?. If m = (left panels), 
both FS states are non-trivial, but with different topological charges {Cps = ±1 or ±2). 
The appearance of Cps=^'^ phase in this case is due to the B^kzk\ term of -B(k). At the 
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boundary between distinct FS states, the fermi spheres should be connected, and the Cps 
becomes ill-defined. The non-trivial FS states may become important for the topological 
superconductivity [27] . 

Even without the exchange splitting, we can expect somethings particular for such novel 
compounds. First of all, we will expect quantum spin Hall effect in z-oriented NasBi thin 
film (or NaaBi/NaaSb quantum well). Due to the quantum size effect, the is further 
quantized, and in general the 2D band structures of NasBi thin film will be fully gapped. 
Then depending on the number of band inversions associated with the subbands, the system 
should crossover between trivial and non-trivial 2D insulators oscillatorily as a function of 
film thickness [28]. Our estimated first critical thickness of NasBi is 35A, below (above) 
which the film is a trivial (non-trivial) insulator. Secondly, we will expect giant diamagnetism 
of 3D massless Dirac fermion [29]. The diamagnetic susceptibility, x(£:)~log^, should diverge 
logarithmically when the chemical potential approaches the MFP (i.e. e~0) [29J, much 
stronger than that in narrow gap semimetals like Bismuth. Thirdly, we will also expect linear 
quantum magnetoresistance (MR) as proposed by A. A. Abrikosov [30j. In conventional 
metals with closed fermi surface, the MR should behavior quadratically at low field and 
saturate at high field. However, for the 3D massless Dirac fermionic gas, the MR will have 
linear field dependence if only the lowest Landau level is occupied. This idea have been 
examined for the Ag2Te both experimentally and theoretically [301 EI]) where the Dirac 
type energy dispersion is not obvious. Having NasBi with MFP, it will be straightforward 
to check the quantum MR proposal [SO]- Finally, it will be very interesting to study the 
possible simultaneous symmetry-breaking introduced by interactions. 

We acknowledge the supports from NSF of China, the 973 program of China, and the 
Hundred Talents Project of the Chinese Academy of Sciences. 

METHODS 

First-principles calculations are preformed based on the plane-wave ultra-soft pseudopo- 
tential method, using the generalized gradient approximation (GGA) for the exchange- 
correlation functional. The calculations based on hybrid functional (HSE) ^I] are further 
supplemented to check the band gap. The cutoff energy for wave-function expansion is 25 
Ry, and the k-point sampling is 12x12x6. The experimental lattice parameters are 
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used, and the calculations are well converged with respect to above settings. The projected 
surface states are obtained from the surface Green's function of semi-infinite system, simi- 
lar to the method used for Bi2Te3 and Bi2Se3 [121 123]. For this purpose, we construct the 
maximally localized Wannier functions (MLWF) from the first-principles calculations. 

The low energy effective hamiltonian is derived from the theory of invariants in a similar 
way as that for Bi2Se3, Bi2Te3 and Sb2Te3 [26]. The first-principles calculations indicate that 
the wave-functions of low energy states at F point are mostly from the Na-3s and Bi-6px,y,z 
orbitals. For the low energy Na-3s at F, about 65% characters are from Na(l)-3s and 35% 
from Na(2)-3s. Since the system has inversion symmetry, the wave-functions at F point have 
definite parity. There are two formular units in each unit cell, and they form bonding and 
anti-bonding configurations. Finally the linear combination of atomic orbitals read: 

\S^) = ^{\Na,s)±\Na,s)), 
\P^) = ^i\Bt,p^)T\Bt\pa)), 

where Na(Bi) and Na'(Bi') are related by inversion symmetry. The superscript ± indicates 
the parity, and a is p^, Py or p^. 

Taking \S^) and {P^) as new atomic orbitals, and including SOC effect, spin and orbital 
angular momenta will be coupled, and the standard eigenstates can be written as |S'f , ±^), 

2 

|P|^,±|), \Pt,±l), \Pf,±l), where the subscript indicates the total angular momentum 

2 2 2 

J. Different from the case with Zinc-blende structure (such as HgTe), here the heavy- hole 
state \Pt, ±1) and light-hole states \Pt, are no longer degenerated at F point, because 

2 2 

Bi atoms are sandwiched by Na(2) atoms along ;z-axis, and Bi pz orbital is lower than p^^y 
orbitals. Furthermore, under the Dqj.^ symmetry, the light-hole state |P|^,±^) and split-off 

2 

state |Pf ,±|) will mix further to form the new eigen states: |P^,±|) and \Pl^,±^) [26] . 
Nevertheless these mixed states are not relevant to our discussions for the MFP. The band 
inversion and their crossings along T-A line can be described by four states |S']^,±^) and 

2 

|P3",±|). Different from Bi2Se3 [22], where all the four basis have the same jJ^H^, here 
they have two different values of | and |, respectively. This difference is essential to the 
existence and stability of 3D Dirac points observed here. 

The parameters in the effective hamiltonian are fitted to the first-principles band struc- 
tures. Their values are: Co=-0.06382 eV, ^1=8.7536 eVA^, C2=-8.4008 eVA^, Mo=-0.08686 
eV, Mi=-10.6424 eVA^, M2=-10.3610 eVA^, A=2.4598 eVA. 
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